PROCEEDINGS 

OF SCIENCE 



Magnetic field production after inflation. 



in 
o 
o 

(N 
Or 

(N 



> 
^f" 

o 

On 
O 

in 
o 

D 



X 



Andres Diaz-Gil 

Departamento de Fisica Tedrica 

Universidad Autdnoma de Madrid, Cantoblanco, 28048 Madrid, Spain 



E-mail: andres . diazgil@uam . es 



J. Garcia-Bellido 

Departamento de Fisica Tedrica 

Universidad Autdnoma de Madrid, Cantoblanco, 28048 Madrid, Spain 
E-mail: 



juan . garciabellidotjjcern . ch 



M. Garcia Perez * 

Institute de Fisica Tedrica UAM/CSIC 

Universidad Autdnoma de Madrid, Cantoblanco, 28048 Madrid, Spain 
E-mail: jnargarita . garcia@uam . es 



A. Gonzalez-Arroyo 

Departamento de Fisica Tedrica 

Universidad Autdnoma de Madrid, Cantoblanco, 28048 Madrid, Spain 



E-mail: antonio . gonzalez-arroyo@uam. es 



We study the electromagnetic field production during preheating after hybrid inflation in a model 
with the field content of the Standard Model, coupled to a singlet inflaton. We find that very 
soon after symmetry breaking our system enters a regime of kinetic turbulence, characterized by 
a self-similar behaviour of the energy spectra and a power-like dependence on time of the inflaton 
and Higgs field variances. 
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1. INTRODUCTION 

There is a large amount of indirect observations of the presence of large scale cosmological 
magnetic fields (LSMF) with scales that go from the size of galaxies to the size of superclusters, 
with a relatively big strength [[I]]. It is believed that these magnetic fields were present in the form 
of a magnetic seed before the epoch of structure formation and that they were amplified to its 
present magnitude. Although there are several attempts to explain these magnetic field seeds, their 
origin is still a mystery It is tempting to speculate whether they could be generated within the 
inflationary paradigm Our present work is along these lines. Taking the end of inflation at the 
electroweak scale, we propose to place the origin of LSMF in the very out of equilibrium epoch 
called the preheating era. Our first aim is to throw some light into the process of electromagnetic 
field production during this period. 

1.1 The model and the method 

We consider a hybrid inflation model which has the field content of the Standard Model (SM) 
gauge-scalar sector plus a singlet scalar field, the inflaton, coupled to the SM Higgs: 



with <J> = i(0o 11 + the Higgs field, % the inflaton and G£ v = - d v A£ + g w e ahc A b ^A c v 

and Fj v = d^By — d v B^ the field strengths of the SU(2) and U(l) gauge fields respectively. The 
covariant derivative is: = — ^g^A^Za — ^g Y B^, with g w the SU(2) gauge coupling and g Y 
the hypercharge U(l) coupling. The scalar potential includes a coupling to a massive inflaton: 

V(O,z)=V o + ^feV-m 2 )|0| 2 + ^| 4 + ^V, (1-2) 

where \(f>\ 2 = 2Tr <£'<!>, \i is the inflaton mass in the false vacuum and m = VXv\ v = 246 GeV. 
During the period of hybrid inflation the inflaton field is dominated by its homogeneous mode %o- 
Inflation ends at t = t c when this mode slow -rolls below % c = m /g where the Higgs becomes mass- 
less. Around this time Xo = Zc(l — Vm(t — t c )), with V the inflaton dimensionless velocity. After 
t c a negative time-dependent mass-squared of the Higgs is induced. During the subsequent evolu- 
tion, the low momentum modes of the Higgs grow exponentially in a process known as "tachyonic 
preheating" [Q]. Preheating in this model has been proposed in Refs. [|], ||] as a mechanism for elec- 
troweak baryogenesis, an appealing possibility that has been numerically investigated in [Q, ^, 

To compute the time evolution of the system after inflation ends, we have used the classical 
approximation (details can be found in ® - see also [§]). For suitable parameter choices, the growth 
of the Higgs field low momentum modes takes place before non-linearities in the Higgs potential 
and the coupling to the gauge fields have become relevant. Thus, one can exactly compute the 
quantum evolution of the Higgs field. It is such that, after a while, the tachyonic modes of the Higgs 
evolve like classical modes. This initial stage is very fast and therefore the remaining degrees of 
freedom can be assumed to remain in their initial quantum vacuum (ground) state. Once sufficient 
number of Higgs modes have become classical we substitute the exact quantum evolution by the 
classical one. For consistency, this substitution should be done before non-linearities have become 
relevant. The advantage of this approximation is that it allows to treat the full non-linear evolution 
of the system by lattice techniques. 
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2. RESULTS 



The results presented here correspond to the particular choice of parameters denoted as Model 
Al in Ref. For this model the inflaton velocity V = 0.024, g 2 = 2A, g w =0.05, and the Higgs to 
W boson mass ratio is given by m H /m w =4.65. Here we also introduce the U(l) hypercharge field, 
fixing g Y to reproduce the experimentally determined Weinberg angle. This choice of parameters, 
although far from the physical value of m w , allows an optimal control of lattice and finite volume 
artefacts for the lattices sizes reachable with our present computer resources [§]. Results for phys- 
ically relevant parameters will be presented elsewhere The chosen values of the lattice cut-off 
are ma= 1.31 (NS=32), ma = 0.87 (NS=48) and ma = 0.655 (NS=64). The time-like lattice spacing 
a t must be smaller than the spatial one for the stability of the discretized equations of motion. We 
have chosen a/a t =40 and 80 and checked stability of the results under this change. 

2.1 Means and Spectra 

In Fig. |l| (Left), the time evolution of (|0|/v) and (x/Xc) illustrates the process of symmetry 
breaking: the Higgs starts in the false vacuum at t = t c and quickly evolves into the symmetry 
broken phase. Fig. [T] (Right) shows the late time behaviour of the different components of the 
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Figure 1: Left: The time evolution of (|0|/v), (x/Xc)- Right: Late time behaviour of the different compo- 
nents of the energy, normalized to the total energy, for NS=48. 



energy, normalized to the total energy, for ma=0.87. The evolution is such that kinetic and potential 
energies virialize as for harmonic oscillators ((kinetic) = (potential)). 

We are interested in the behaviour of electromagnetic fields long after symmetry breaking. It 
is not obvious that our approximations remain under control for such late times. The validity of the 
classical approximation rests on the fact that the tachyonic infrared modes dominate the dynamics. 
This cannot happen if ultraviolet modes are strongly excited when the physically relevant processes 
take place. A way to expose the range of relevant momenta is to analyze the Fourier spectra of the 
(gauge invariant) energy densities. Fig. § shows, as an illustration, the time evolution of the spectra 
of the inflaton+Higgs contribution to the kinetic energy and of the electric component of the SU(2) 
gauge field energy. To enhance the effect of lattice artefacts for large momenta, we plot k 2 \E(k)\. 
Note that the presence of the lattice cut-off truncates the range of allowed momenta. The number, 
N(k), of lattice momenta of norm k should be proportional to k 2 in the range of relevant momenta 
for the cut-off not to distort the physical spectrum. Deviations of N(k) from the k 2 behaviour appear 
for k/m > 2.5, 3.5, 4.5 for NS= 32, 48, 64 respectively. It is clear from the figure that the NS=32 
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lattice is off. NS=48 starts showing deviations, for physically relevant momenta, for times above 
mt ~ 100. These deviations, however, seem to remain small up to much larger times. Full checks 
of cut-off independence by going to larger lattices will be done in the near future. Nevertheless, 
given these results, our approach to study the evolution up to relatively large times, mt ~250, seems 
viable. The time ranges in this study are four times larger than those in Ref. [^]. 




Figure 2: Spectra of the inflaton+Higgs kinetic energy and SU(2) electric energy. 



2.2 Electromagnetic field production 

In order to analyze the production of electromagnetic fields we have, first, to extract the U(l) em 
content of the SU(2)xU(l) fields in the Lagrangian. Details of the lattice implementation will be 
given elsewhere In terms of the SU(2) link, U^(x), and the hypercharge link, B^(x), one can 
express the field associated to the Z boson as: 

with = a, ciq = a t and gz the Z boson coupling. The field strength of the U(l) em field is 

obtained from: F^{nf = 1 - cos {jr^rF^ v {n) - f/ v (n)) , with F^ v (n) = Z ll (n) -Z v (n + fi) - 

Z^{n + v) — Z v (n) and F^ v (n)=B ii {n)—B v {n + jjL) — B ii (n + v) — B v {n). In the continuum limit 
— > a^a v eF^, with F™ the continuum electromagnetic field strength. 
In computing the field strength of the Z boson, F^ v , from Eq. ( |Q| ) there is a complication that 
arises due to the possible discontinuous behaviour of <&' (n)/\(j>(n)\ at points where the Higgs field 
is zero. At such points the Z boson and the electromagnetic field strengths can show a (lattice) delta- 
function behaviour. This gives rise to the appearance of non-scaling spikes in the corresponding 
energy densities which affect also spatially averaged quantities. Although configuration averaging 
smears out the effect of the spikes, they give rise to a very noisy signal. One way to improve 
the signal is to eliminate from the spatial averages the points where the Higgs norm gets closer to 
zero than a certain value \<p c \. Taking \(j> c \ = 0.08 considerably reduces the noise and gives results 
compatible within errors with those obtained from larger statistics sets. The appearance of these 
spikes is rather common in the first stages of the evolution where the Higgs field gets rather often 
close to zero. It becomes more rare as time evolves, although occasionally they can also occur 
at times at large as mt = 200. It is atttactive to speculate that these zeroes could be associated to 
sphaleron-like configurations, a point that deserves further investigation [^]. 
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Figure 3: Left: Late time behaviour of the electric and magnetic components of the U(l) e m energy, normal- 
ized to the total energy. Right: Longitudinal and transverse parts of the electric component. 



Figure |3| shows our final result for the late time behaviour of the electric and magnetic com- 
ponents of the U(l) em energy, normalized to the total energy. Results are shown for two different 
lattice spacings, agreeing well within errors. A detailed analysis of the properties of the electromag- 
netic field will be deferred to a future paper [Bp. Here we only want to point out that a separation 
between the longitudinal and transverse components of the electric field can be easily performed 
by going to Fourier space. The result is shown in Fig. |5| (Right). The longitudinal component 
decreases with time. This is the expected behaviour at large times since the system evolves into a 
plasma with charge density going to zero. 
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Figure 4: Left: The late time behaviour of the variance of the Higgs and infiaton fields. Right: Spectra of 
the SU(2) electric part of the energy showing a self-similar scaling according to Eq. (2.2). 



2.3 Kinetic turbulence 



According to Ref. Q10|], systems like the one we study evolve towards thermalization via a 
period of kinetic turbulence. This regime is characterized by a specific behaviour of the variance 
of the scalar fields together with a particular scaling of the energy spectra. To be concrete, one 
expects for the scalar fields: ((f) 2 ) — (0) 2 ~ t~ y with v = 2/(2m— 1) in three spatial dimensions 
and for m-particle interactions. At the same time the energy spectra show a self-similar behaviour: 



n(k,t) = t q no(kt p ), 



with n(k,t) 



E(k,t) 



(2.2) 



where q = 3.5p and p = I /(2m — 1). Ref. [ |10| ] analyzes models containing only scalar fields. 
The characteristic times for the kinetic turbulence regime are mt > 3000 after inflation ends. We 
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also find that our data enter a regime of kinetic turbulence, but this happens at much earlier times 
mt > 50, presumably due to the larger number of degrees of freedom. The variances of the Higgs 
and inflaton fields are shown in Fig. |](Left). They follow the expected behaviour with m =4 for the 
Higgs and m = 3 for the inflaton. Fig. |] (Right) shows the spectra of the electric part of the SU(2) 
energy for several times from mt = 55 up to mt = 250. The data remarkably follow the self-similar 



behaviour of Eq. (2.2) with p = 1.1/7, very close to the value expected for m = 4. The early onset 
of kinetic turbulence is good news since it allows extrapolation of the time evolution beyond the 
limitations set by our numerical approach. 

3. CONCLUSIONS 

Making use of the classical approximation, we have studied the production of electromagnetic 
fields during preheating after a period of low-scale hybrid inflation. We present here a feasibility 
study showing how our approach reliably describes the evolution up to rather large times, m?~250. 

We have extracted the total electric and magnetic t/(l) em energies. Concerning electric fields 
we observe a damping of the longitudinal component of the field strength, corresponding to the 
gradual neutralization of charged particles in the primordial plasma. The total magnetic field adds 
up the contribution of radiation and the primordial magnetic seed. The separation of the latter 
will be deferred to a future publication [Q]. It is particularly interesting the fact that, very soon after 
inflation ends (mt > 50), the SU(2) energy spectra show a self-similar behaviour (see Fig. (Q)). This 
is a signature of kinetic turbulence. We still have to check that this seed gives rise to a magnetic 
field correlated over large scales (comparable to the horizon), and that it remains correlated for a 
sufficiently long time. 
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